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O ■ Abstract 

, We discuss simulations with different lattice Dirac operators for sea and valence quarks. A 

ry^ | goal of such a "mixed" action approach is to probe deeper the chiral regime of QCD by enabling 

■ simulations with light valence quarks. This is achieved by using chiral fermions as valence quarks 

. while computationally inexpensive fermions are used in the sea sector. Specifically, we consider 

\ Wilson sea quarks and Ginsparg- Wilson valence quarks. The local Symanzik action for this mixed 
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theory is derived to 0(a), and the appropriate low energy chiral effective Lagrangian is constructed, 

CN ■ including the leading 0(a) contributions. Using this Lagrangian one can calculate expressions for 
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physical observables and determine the Gasser-Leutwyler coefficients by fitting them to the lattice 
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I. INTRODUCTION 



In order to extract predictions of QCD from numerical methods with controlled systematic 
errors, a lattice formulation is required for which the sources of deviations from QCD are 
understood and are under control. A significant source of systematic errors for present day 
lattice simulations are the light quark masses. Even the most powerful computers today do 
not allow simulations with up- and down-quark masses as light as realized in nature. Instead 
one simulates with heavier quark masses and fits the analytic predictions obtained in chiral 
perturbation theory (xPT) to the data. The free parameters in the fit are the low energy 
constants ofyPT |l[ , and once they are determined an extrapolation to small quark masses 
is possible |2j, |3( . Still, to perform the chiral extrapolation the quark masses must be small 
enough so that is applicable. In practice one would require that next-to-leading order 
(NLO) xPT describe the data reasonably well. 1 
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The data do 



The present lattice simulations do not meet this requirement 
not show the characteristic curvature predicted by NLO x?T. In fact the data show a 
rather linear behavior which either means that higher orders in the chiral expansion are not 
negligible or worse, one is not in the chiral regime at all (see Claude Bernard's part in j^J]). 
In either case, simulations with lighter quark masses are required in order to apply %PT 
with confidence. 

Lattice simulations with light fermions, especially sea quarks, are computationally de- 
manding and the numerical cost increases substantially with decreasing quark masses. Re- 
alistically only the least expensive fermions, Wilson and Kogut-Susskind, can be used on 
sufficiently large and fine lattices. Lattice fermions with better chiral properties are still too 
expensive to be used as sea quarks, and this situation is not likely to change in the near 
future. It is nevertheless expected that the next generation of TFLOP machines will make 
it possible to generate a few sets of unquenched configurations with sea quarks light enough 
to be in the chiral regime. 

To obtain more information from these configurations they should (and will) be analyzed 
with various different valence quark masses, i.e. by studying PQ QCD. By including lattice 
measurements with lighter valence quarks it is possible to penetrate further the chiral regime 



1 At next-to-next-to-leading order (NNLO) many new unknown parameters enter the chiral Lagrangian, 
which greatly reduces the predictive power of xPT- 
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of QCD. This leads to more data points and would allow more reliable fits of PQ xPtQ 
to the lattice data 2]. The reach of such simulations, however, is limited. The cost of light 
valence quarks also increases with the decreasing mass, and can become prohibitively high 
for quark masses that are still not very small. This is particularly true for Wilson fermions 
because of the explicit chiral symmetry breaking by the Wilson term. 

An interesting idea for probing the chiral regime is to use different lattice fermions for 
the valence and sea quarks. In particular, by choosing lattice fermions with good chiral 
properties for the valence quarks, the valence quark mass can be made much smaller than 
in ordinary PQ simulations. A central goal of this strategy is the same as of PQ QCD 
- to explore a larger portion of the chiral regime by extracting more data points from a 
given set of unquenched configurations (see fig. [TJ. This should result in more reliable 
estimates for the low-energy constants of %PT at NLO, the Gasser-Leutwyler coefficients. 
Furthermore, one might expect to reduce the size of explicit chiral symmetry breaking by 
using Ginsparg- Wilson fermions at least for the valence quarks. This is a computationally 
affordable compromise of the lattice theorist's ideal of using Ginsparg- Wilson fermions for 
both valence and sea quarks. 

In this paper we construct the low-energy chiral effective theory for a "mixed" lattice 
action, with explicit dependence on powers of the lattice spacing a by first constructing 
the appropriate local Symanzik action. There are several reasons for taking this approach. 
First, the defining non-orthodox feature of the mixed action approach - the use of differ- 
ent Dirac operators for the sea and valence sectors - is purely a lattice artifact. This is a 
consequence of the fact that by construction all proper lattice fermions reproduce the same 
continuum physics, and therefore all mixed lattice theories reduce to PQ QCD in the con- 
tinuum limit. An expansion in a is thus a natural tool to investigate potential peculiarities 
of the mixed action formulation. Second, a theoretical understanding of the a-dependence 
in lattice simulations can guide the continuum limit, or allow the extraction of physical 
information directly from the lattice data, without taking the continuum limit first. Third, 
%PT provides a useful framework for studying the chiral symmetry breaking due to the dis- 
crete space-time lattice. Effective theories of this type have been studied in several similar 



contexts 



ay 



hi. 



What is dubbed here "mixed action methods" refers to a class of lattice theories corre- 
sponding to different choices of Dirac operators for the valence and sea quarks. In the next 
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FIG. 1: Qualitative representation of the space of quark masses. The "chiral regime", where PQ 
XPT can be applied, is the quarter- circular region. The upper right rectangle, limited by the dashed 
line, describes the part of the space covered by present simulations. As current data suggest, there is 
very little or no overlap between that rectangle and the chiral regime. It is expected that improvement 
in algorithms and computer power will allow reducing the sea and valence quark masses in PQ 
simulations, as is represented by the enlargement of the previous rectangle. It is possible that 
the chiral region will be penetrated by such simulations, as shown by the small section of overlap 
between the enlarged rectangle and the chiral region. Finally, using chiral valence fermions in a 
mixed action simulation would make it possible to extend the reach of simulations significantly in 
the direction of lighter valence quarks. 

section we use a fairly simple example to illustrate the general framework of mixed lattice 
theories. We consider Wilson fermions for the sea quarks, together with valence fermions 
that satisfy the Ginsparg- Wilson relation. To describe the lattice action close to the contin- 
uum limit we construct the local Symanzik effective action up to 0(a). The usual arguments 
used in the formulation of are then applied to this effective action. This leads to a 
chiral expansion in which the dependence on the lattice spacing is explicit. 
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II. THE CHIRAL EFFECTIVE ACTION 



A. Lattice action 

In the following we always consider a hyper-cubic Euclidean space-time lattice with lattice 
spacing a. We assume either an infinite lattice or a finite lattice large enough that one can 
safely ignore finite volume effects. 

The mixed lattice action that describes Nf Wilson sea and Ny Ginsparg- Wilson valence 
fermions has the structure 



>M 



Sym[U] + S w [ip s ,ip s ,U] + S GW {ip v ,ip v , U]. (1) 



U denotes the gauge field defined on the links of the lattice, ips C0s) are the sea quark 
(anti-quark) fields and ipy (ipy) denote vectors with Ny anti-commuting valence quarks 
(anti-quarks) and Ny c-number- valued ghost quarks (anti-quarks). 

The precise choice for the gauge field action Sym is irrelevant in the following, so we leave 
it unspecified. For the sea quarks we choose the Wilson action 12j, given by 



S w = a^^ipsiDw + m Sea )^s(x) , (2) 



X 



D w = \ K( V* + V„) - a r V* V„) } , (3) 

where mg ea denotes the Nf x Nf quark mass matrix in the sea sector and r the Wilson pa- 
rameter. V* , V M are the usual covariant, nearest neighbor backward and forward difference 
operators. 

The action for the valence and ghost quarks is given by 

Sow = a 4 ^2?py {D GW + m Va i (l - \aD GW )} ipy(x) . (4) 

X 

The valence and ghost quark masses are contained in the 2Ny x 2Ny mass matrix myai of 
the form m Va i = diag(M Va i, M Va i) where M Va i is an N v x N v matrix (i.e. each valence quark 
has a corresponding ghost field with the same mass). The Dirac operator Dqw is assumed 
to be a local operator satisfying the Ginsparg- Wilson relation [3| 

lbD G w + £W75 = aD G wj 5 D Gw . (5) 



Both the fixed-point Dirac operator and the overlap operator satisfy this relation 14j,[l5|, 16^ . 



For the following discussion, however, there is no need to specify D G w an Y further. 



B. Flavor symmetry of the lattice action 



When mg ea = 0, myai = 0, and r = 0, the flavor symmetry group of Sm is 

SU(N f ) L ® SU (Nf)n ® ^[/(iVyliVv)^ <g> 5£/(iVv|JVv)fl. (6) 

To see this it is convenient to write Saw an d SV i n terms of chiral components. The right- 
and left-handed sea quark fields are defined with the usual projectors |(1 ± 75). For the 
valence and ghost fields, one first defines the hermitian operator 

75 = 75(1 - aD GW ) , (7) 

which is unitary as a consequence of ( |5|) . Valence right- and left- handed fields are now 
defined by H 

Vv,fl = t/Vt^-Ts), ipv,R = ^(i + l5)ipv, (8) 
Vv,l = V'y^l + Ts), Vv,z = -ji^-l^v- (9) 
The fermionic actions can now be rewritten as 

S w = a±Y/$s, L \l» (V, + V;) *P s>L (x) + ^§7^ (V, + vy ^(x) (10) 

a; 

+ ^S,L (™Sea - §<""V M V*) ^ S ,r( X ) + ^S,R ( m L ~ \ ar ^ M V ^) 1>S,l( X )> 

and 

Sew = a 4 ^ V)L Vv,l(x) + £>gw Vv,i?(x) 

x 

Here, for reasons that will become clearer shortly, we consider ms ea , my a \ and r to be matrices 
in flavor space, and identify the parameters that appear between right-handed anti-quarks 
and left-handed quarks as their hermitian conjugates. 

Clearly, when ms ea = "7,Vai = r = 0, eqs. (jlOllllj) are invariant under independent global 
rotations of the left-handed and right-handed components of all quark fields: 

*Px, x -> gx^xa, ^x, x ~> ^x,x9x,x> X = S,V, x = L,R, (12) 
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where gs,L and g$,R are in SU(Nf), while gv,L and gy.n are in S77(AV|iVV). We conclude 
that flavor transformations belonging to the group in eq. © are symmetries of the action 
eq. (JTJ broken by the mass terms my a i and msea and the Wilson term r. 

It is convenient to treat these symmetry breaking parameters as spurion fields, i.e. as- 
suming the transformation rules 

m Yal -> g V ,L ™Val 9v,R , m l/ a l flV.fl m Val #V,L > 

™Sea -> C/5,L ™Sea > "4a ~» "4ca g\ L , (13) 

r -> 9s,Lrgl iR , r f -> g s ,R^ g\ L . 

The mixed action 5m, even with non- vanishing mass and Wilson terms, is invariant under 
the combined transformations eqs. (jl2H13|) . 

To complete this part we note that for a transformation to be a symmetry of the theory it 
must also leave unchanged the integration measure in the functional integral. It is a simple 
matter to show that the measure for the sea Wilson fermions is invariant under the global 
rotations in eq. (|T2|). The situation for the valence quarks is not quite so simple because 
of the operator 75 in the chiral variation. It turns out, however, as has been shown in ref. 
[l8|. that the measure is indeed invariant under the symmetry transformations considered 
here - the flavor non-singlet transformations. The last statement can be extended to the 
full valence sector, including the ghost fields. 2 



C. Symanzik action 

We construct Symanzik's local effective theory which, close to the continuum, describes 
the same long-range physics as the discrete lattice action well below the momentum cutoff 



1/a |21|,|22,|23|. 



Since the continuum action 5*5 is designed to reproduce the same long-range correlation 
functions as the discrete lattice action Sm, it must have the same symmetries [eq. (jHJ)] as 
the underlying theory. Up to 0(a), the quark operators that enter are of mass dimensions 



2 The symmet ry g roup that we write here is not the true symmetry group of the quantized theory. As 
discussed in [id Epf. the presence of ghost fields in the functional integral leads to constraints on the 
allowed symmetry transformations. However, for the derivation of the correct chiral Lagrangian it is 
possible to use the symmetry group in eq. JBJ [iof . 
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3,4, and 5, which include only quark bilinears. Moreover, the independent symmetry trans- 
formations acting separately on the sea and valence sectors requires that the quark bilinears 
do not mix the sectors. This implies that up to 0(a), the fermionic operators in S$ (as in 
Sm) are of two types - one built of sea quarks only and one of valence quarks. 

It is straightforward to write down the 0(a) Symanzik action Ss using previous results 



concerning Wilson fermions |25l. l26| and Ginsparg- Wilson fermions [17 1. The details of the 
analysis are deferred to appendix [X] ~~ here we only quote the result (for the fermionic part 
of the action): 

S s = J d 4 x [ip v (D + mvaOVv + i>s( D + m Sea )ips + ac S w^s a ^ F ^s] + 0{a 2 ). (14) 

mvai and ms ea are renormalized masses. Two consequences of the exact chiral symmetry 
of the massless Ginsparg- Wilson fermions are: (a) there is no Pauli term ipa^F^ijj for the 
valence sector, and (b) the valence quark mass is only multiplicatively renormalized. No 
symmetry protects the Wilson sea quarks from getting an additive correction of the order 
of the cutoff I /a. 

It is useful at this point to collect the quark fields in a single quark field vector and 
rewrite eq. (fT4^) as 

S s = J^(D + m )^+^ac sw a^F^ + 0(a 2 ) , (15) 

( i's\ (rh Se& \ (ac sw \ 

* I I ■ m = 1 . r/c S w =1 I • (16) 

D. Symmetries of S$ and xPT for the mixed action 

We now turn to the construction of a low-energy effective theory for the "underlying" 
Symanzik action in eq. (|15j). The method is completely analogous to the construction of 
the chiral Lagrangian in QCD j^]. 3 The idea is that the spontaneous breaking of the 













, m = 




, ac S w 






^ ^Val J 





3 It should be noted that as of yet the construction of PQ xPT from PQ QCD is not as well justified as 
the standard derivation of x?T from QCD. The arguments of the latter cannot be trivially extended to 
PQ QCD because they rely, in part, on the existence of a Hilbert space of physical states with a positive 
definite norm, which is absent in the presence of ghost fields. The same is also true for the mixed action. 
The validity of the yPT for the mixed action is thus on the same footing as that of PQ xPT, which has 
been discussed in 201. 
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approximate chiral symmetry gives rise to light pseudo Goldstone bosons, the light mesons, 
which at low energies are the only relevant degrees of freedom. The effective action is written 
in terms of local interactions of the pseudo Goldstone fields. Since the pseudo Goldstone 
bosons interact weakly at low energies, the action can be organized in a perturbative series. 
All observables calculated are expanded in two small parameters, 

X X X 

where p is the light meson momentum, A x ~ lGeV is the chiral symmetry breaking scale, 
and m and a stand for matrix elements of 

m = 2B m, a = 2W ac sw . (18) 

B and Wq are dimensionful low-energy constants that appear in the effective theory at 
leading order [see eq. (|B1)) ]. They depend only on the high-energy scale A x , and dimensional 
analysis reveals that the perturbative expansion is in fact in m q /A x and aA x . We follow 
ref. |9[ and choose the {e, 5} terms as leading order (LO) and the {e 2 ,e<5} terms as NLO in 
the effective Lagrangian, dropping 0(a 2 ) contributions. The underlying hierarchy consistent 
with this ordering is {e, 5} ^> {e 2 ,e5} ^> 5 2 , and the last inequality also implies e 5. 
This ordering is chosen for convenience and is somewhat arbitrary. In practice, the double 
expansion should be organized according to the actual relative sizes of the quark masses and 
the lattice spacing. 

The effective Lagrangian is constructed from all operators that respect the symmetries 
of the underlying action S$- The compact notation of eq. (|15|) makes it easy to see that the 
symmetry group of Ss to 0(a) is 

SU(Nf + N V \N V ) L ® SU{N f + N V \N V ) R . (19) 

This symmetry group (treating m and a as spurion fields) is the same as that of PQ 
QCD with Wilson fermions. Indeed the effective action Ss in eq. (|15j) is the same as the 
effective Symanzik action for the PQ QCD Wilson action of ref. |9j, with a specific choice of 
acsw that has support only in the sea sector. This fortunate similarity between the mixed 



4 Though the notation might obscure this fact, m and a both have mass dimension 2 - they are the leading 
contributions to the squared mass of the pseudo Goldstone boson. We nevertheless use this notation as 
it makes the dependence on the quark masses and the lattice spacing more transparent. 
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action theory and the PQ QCD Wilson theory implies that the mixed low-energy chiral 
effective action has the same structure, at 0(a), as the action of Wilson x?T ; introduced 
in ref. \f^, with the restriction that acsw vanishes in the valence-ghost sector. The chiral 
Lagrangian to 0(a), which describes both W%PT and the mixed effective theory, is provided 
in appendix FBI 



E. Application: Meson mass 

In this subsection we give an example for the use of the mixed chiral Lagrangian. For 
simplicity we take all the sea quarks and all the valence quarks to be (separately) degenerate, 
and the Wilson parameter to be a flavor singlet in the sea sector. This amounts to setting 

m = diag (m So a, wval) , a = diag (a, 0) . (20) 

The number of sea quark flavors is taken to be Nf = 3. We consider the expression for the 
mass of the flavor charged meson with valence quark flavor indices AB (A ^ B) to NLO. 
Using the relation between the mixed chiral effective theory and W^PT, one can obtain the 
result straightforwardly by taking the mass formula from ref. j^] with the values for m and 
a given by eq. (J20J) . We find 



( M AB) NLO = ^Val + 4 gyS^2 [™ Val ~ ™ Sea ~ " + ( 2 ™Val ~ ™Sea ~ o) In (m V al)] (21) 

[(L B " 2L 8 )m V ai + 3(L 4 - 2L 6 )m Sca + 3(W 4 - W 6 )a] . 



Bl 



r- 

Here, the parameters Lj are the Gasser-Leutwyler coefficients, and W4 and Wq are additional 
low-energy constants that enter the chiral Lagrangian at NLO. Note that for a = the 

28, 291) is recovered. 



expression for PQ x?T (calculated in 

Eq. (12 lj) demonstrates the analytic connection between QCD and the simulated mixed 
action theory. It shows the latter to be a calculation with controlled systematic errors. From 
fitting the equation to the appropriate data from numerical simulations one can obtain an 
estimate for the linear combinations of Gasser-Leutwyler coefficients that appear in it. 

Examining eq. (|21j) one can also appreciate the potential advantage of using a mixed 
lattice action. In simulations using Wilson fermions in both sea and valence sectors, an 
equation similar to eq. (}2~T]) holds (see eq. (|B5(l ). In that case, the range of valence quark 
masses that can be simulated might be too small to convincingly show the curvature coming 
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from the quadratic dependence and the logarithms that enter at NLO. By using Ginsparg- 
Wilson fermions for the valence quarks one can vary the valence quark masses over a wider 
range. The expected NLO curvature, on which the extraction of the Gasser-Leutwyler 
coefficients depends, is consequently much more likely to be seen. 

Finally, comparison with the result for W^PT, eq. (jB5|) . reveals that the latter depends 
on twice the number of Wi coefficients. This is fortunate for the mixed theory as it makes 
the predictions of the effective theory less dependent on parameters that have no particular 
relevance to QCD. 

To understand this simplification in the expression for the meson mass, consider the 
relation between the symmetries of the mixed action and those of PQ Wilson action. On the 
one hand the massless mixed theory has exact chiral symmetry in the valence sector, which 
the Wilson action does not. On the other hand, the valence and sea quarks of the Wilson 
action have the same type of Dirac operator which allows mixing between the sectors - a 
transformation which is not a symmetry of the mixed action formulation. At 0(a), however, 
the breaking of the sea-valence symmetry in the mixed theory does not yet show up, and 
thus the simpler expressions arise due to the larger chiral symmetry. 

III. SUMMARY 

In this paper we discuss lattice simulations with different fermions for sea and valence 
quarks. As a particular example we have studied here the case with Wilson sea quarks 
and Ginsparg- Wilson valence quarks. Using Symanzik's effective action for lattice theories 
as an intermediate step, we have derived the form of the low-energy chiral Lagrangian for 
the mixed theory to 0(a). The construction shows that simulations with the mixed action 
provide as controlled an approximation to QCD as partially quenched simulations. This is 
to be expected since the mixed action reduces to PQ QCD in the continuum limit. 

The goal of the mixed action approach is similar to that of PQ QCD. The use of smaller 
valence quark masses allows one to probe deeper the chiral regime of QCD and obtain 
additional information on the low-energy constants, the Gasser-Leutwyler coefficients. Fur- 
thermore, the use of chiral lattice fermions in the valence sector, instead of Wilson fermions, 
makes it possible to simulate much lighter valence quarks. This leads to more data points 
obtained on the lattice and consequently to more reliable fits of %PT to the data. 
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Here we have demonstrated the mixed action approach for Wilson sea quarks and 
Ginsparg- Wilson valence quarks taking into account the leading 0(a) contributions. An 
important extension of this analysis is the inclusion of 0(a 2 ) effects. First of all, the lat- 
tice spacing determined by the unquenched configurations is possibly not small enough to 
safely neglect the 0(a 2 ) corrections. If one wants to fit the lattice data directly to equations 
like eq. ()21|) without taking the continuum limit first, the 0(a 2 ) corrections should be in- 
cluded to obtain better fits. Second, the 0(a) effects are generated here only by the Wilson 
sea quarks. Many unquenched simulations are in fact performed with non-perturbatively 
(9(a)-improved Wilson fermions for the sea quarks. The leading corrections for these simu- 
lations are of 0(a 2 ) and need to be computed in order to know how the continuum limit is 
approached. 

While valuable, it should also be noted that the inclusion of 0(a 2 ) effects in the chiral 
Lagrangian framework is likely to be a hard task. The main difficulty arises from the 
many new operators that enter the Symanzik action at this order. Some of these operators 
break Lorentz invariance, while several others break the chiral symmetry and require the 
introduction of additional spurion fields. 

The approach proposed here should be also studied with other combinations for the lattice 
fermions. In particular, the case with staggered sea quarks is interesting, since staggered 
fermions are computationally cheaper. At present the lightest dynamical quark masses are 
achieved with staggered fermions. It is well-known that applying staggered fermions to QCD 
involves a theoretical uncertainty and is possibly uncontrolled. Consequently, predictions 
from chiral perturbation theory for staggered fermions would also serve as a test of this 
discretization method Q. 

Finally, the cost of simulations of a mixed action is roughly the sum of the cost of 
generating a set of unquenched gauge field configurations plus that of analyzing quenched 
simulations with Ginsparg- Wilson fermions. Thus, we can expect that in the near future 
simulations with a mixed action will become feasible. 
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APPENDIX A: SYMANZIK ACTION FOR THE WILSON AND GINSPARG- 
WILSON ACTIONS 

In this appendix we derive eq. (J 14)) for the Symanzik action describing the mixed theory 
to 0(a). As has been stated in the text, to this order the Symanzik action is simply the 
sum of the local effective actions for the valence and the sea sectors. 

The local Symanzik action for Wilson fermions has been derived in (25! [26^ . One first 
lists all the operators of mass dimension no greater than 5, which respect the symmetries of 
the Wilson lattice action (the appropriate power of a is inserted to complete the dimensions 
of terms in the Lagrangian to 4). The operators of dimension 4 (which are a independent) 
make up, by construction, the continuum action of QCD. 

Because the Wilson term explicitly breaks the chiral symmetry, it is expected that the 
quark mass be additively renormalized, and the size of the correction should be of the order 
of the cutoff scale 1/a. Indeed, the only dimension 3 operator is ipip, which appears in the 
action with a coefficient proportional to 1/a and has precisely this effect. 

There are several operators of mass dimension 5 that are allowed by the symmetries. 
Some of these operators can be eliminated using the leading order equations of motion. 
Others have the same structure as the mass and kinetic operators that already appear in 
the QCD action, and have the effect of renormalizing the quark masses and the gauge 
coupling. Finally, a single term is left - the Pauli term: ipcr^F^il). Note that the Pauli 
term breaks the chiral symmetry, and is therefore allowed only because of the Wilson term. 
Putting it all together, the Symanzik action for the Wilson sea sector is 

J d 4 x \$ S {D + m Sea )V>s + ac sw i's a ^ F ^i's] + 0(a 2 ) . (Al) 

where msea is the renormalized sea quark mass, and csw is an unknown coefficient. 

The analysis for the Ginsparg- Wilson valence quarks is similar. This may seem confusing 
due to the fact that some of the chiral projectors on the lattice are written in terms of 75, 
and not 75 as in the continuum theory. However, it has been shown in ref . 0, E3| , that the 
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chiral symmetry of the Ginsparg- Wilson lattice action leads to exactly the same chiral Ward 
identities which appear in the continuum. Hence, by imposing the usual chiral symmetry 
on the Symanzik action, the effective theory correctly reproduces the consequences of the 
lattice chiral symmetry. 

Due to the exact chiral symmetry the valence quark mass gets renormalized only multi- 
plicatively and the Pauli term is absent. Consequently, after considering the renormalizations 
of gauge coupling and quark masses, the Symanzik action for the valence Ginsparg- Wilson 
quarks contains no 0(a) part (the Ginsparg- Wilson lattice action is automatically 0(a) 
improved 



Eq. ((HJ) is the sum of eq. J53J and eq. (jX2j). 
APPENDIX B: W X PT RESULTS 

We present the W%PT Lagrangian which also describes the mixed theory to 0(a). In- 
terested readers should consult ref. |i| for further details on W%PT. We also provide the 
expression for the mass of a flavor charged meson for comparison with the mixed theory 
result. 

The W^PT Lagrangian is constructed out of operators that respect all the symmetries 
of the underlying theory in eq. (|15p. with explicit flavor axial symmetry breaking terms 
constructed out of in and a. As described in the text, the LO Lagrangian is linear in e and 




(A2) 



8: 




(Bl) 



Here the angled brackets stand for the super-trace over the flavor indices: 



(r) = str(r) = J^riiTii, rji 



1 i is a quark flavor index 
- 1 % is a ghost flavor index 



(B2) 



and £ = exp (2iU/f) is a non-linear representation of the meson fields. 
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The NLO Lagrangian is: 

£ 4 = L^dZd^) 2 + L 2 (d^d u ^) (d^d u ^) + L 3 ({dZd^) 2 ) 
+L 4 (dEdE 1 ) (m f E + E f m) + W 4 (dEdE 1 ) (a f E + E f a) 
+L 5 (9E9E^ (m^E + E+m)) + W 5 (dEdE^E + E f a)) 
+L 6 (m f E + E f m) 2 + W 6 (m f E + E f m) (a f E + E f a) 
+L 7 (m f E - E f m) 2 + W 7 (m f E - E f m) <a f E - E f a) 
+L 8 (t^Et^E + YtmYtm) + (^Err^E + EWm) . (B3) 

These Lagrangians describe both the mixed and the PQ Wilson lattice actions. In the mixed 
theory a has support only in the sea-sea sector. We comment that eqs. (|B1|) . (|B3|) contains 
ordinary \PT. Moreover, since the low-energy constants Lj's and W{S are independent of 
m and a and this theory becomes the familiar xPT in the sea-sea sector when a — > 0, the 
Li's are exactly the Gasser-Leutwyler coefficiants of ordinary %PT . 

Next, we provide the W;\PT expression for the mass of the flavor charged meson defined 
in Secion ITTeI We consider the case where 

m = diag (m Sea , ™Vai) , a = diag (a Sea , a Va j) (B4) 

(compare with eq. (J2DI))- O ne obtains: 

/A/T2 \ /~ , ~ \ , ( m Val + aval) r/- , - \ , ^ \ 

( M AbJ NLO =( m Val + CtVal) H 43^2 i( m Val + V al) ~ (^Sea + ^Sea) 

+ (2(m Va i + aval) - (m Sca + a Sea .)) In (m Va i + Oval)] 

8m Va i 



[(L 5 - 2L 8 )m Va i + 3(L 4 - 2L 6 )m Se a + 3(W 4 - W 6 )as«J 
, [(L 5 + W 5 - 2PV 8 )mvai + 3(L 4 - W & )m^} (B5) 



/ 2 

taval 



/ 

To obtain the expression appropriate for common lattice simulations, in which the Wilson 
term is the same for all flavors, one sets av a i = «sea in the last equation. The meson mass 
for the mixed action [eq. (|2~Tj) ] can be obtained by setting ay a i = 0. 
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